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^ I. INTRODUCTION 

The study of exactly solvable models of strongly correlated electrons is important for understanding fundamental 
, aspects of statistical mechanics. It is relevant to many realistic physical systems such as high- Tc superconductors. 
The quantum inverse scattering method (QISM) has been applied to solve various strongly correlated electron models 
0. Two models vifhich have attracted a great deal of attention are the Hubbard model, which is derived from an 
_R- matrix with non-additive spectral parameter, and which was investigated by Shastry 0, and its strong coupling 
\ limit the t — J model. Essler and Korepin established the integrability of the t — J model by obtaining an infinite set 
■ of conserved quantities, and studied the model in the framework of the graded QISM Jp[. A integrable anisotropic 
I ' (g-deformed) supersymmetric t — J model was proposed by Foerster and Karowski [QT Woynarovich ||^ applied 
''O the finite-size correction of the ground state energy to obtain the low-lying gapless excitation spectrum around the 
Ch ground state. Frahm and Korepin obtained the critical exponents for various correlation functions. Furthermore, 
^ the Bethe ansatz solution and conformal properties of the g-deformed t~ J model were studyed by Bariev, Kliimper, 
I I Schadschneider and Zittartz [Q. 

\ By considering other representations of quantum superalgebras, new integrable, strongly correlated electron models 
^""^ ■ of interest, such as the supersymmetric Essler-Korepin-Schoutens (EKS) extended Hubbard model [||, have been 
1 ^ proposed. This gl{2\2) supersymmetric model contains the supersymmetric t — J model as a submodel and can be 
interpreted as the Hubbard model plus moderate nearest-neighbour interactions. The complete solution of the EKS 
model by the algebraic Bethe ansatz has been obtained |^ . The mathematical issue of the completeness of the solution 
l/-^ , has been settled @ , and the physics content of the solution, low lying excitations in particular, has been studied . 

Another non-standard integrable model of strongly correlated electrons is the Bariev chain ||T^ , which has an _R- matrix 
T-H with non-additive spectral parameter as investigated by Zhou |p^ . A generalization of the Hubbard model with Lie 
superalgebra gl(2\l) supersymmetry, the supersymmetric U model, was discovered in and has been investigated 
by several groups |15|-|lq|. An extension of this model, a g-deformed version with quantum superalgebra Uq{gl{2\l)) 
supersymmetry, was also proposed p9| , |20| ]. Thermodynamic properties of the EKS model and the supersymmetric 
U model have been studied using Wiener-Hopf techniques and the critical exponents of correlation functions by 
using methods of conformal field theory [|9| . Recently, the eight-state supersymmetric U model of strongly correlated 
fH ' electrons with the Lie superalgebra gl(3\l) symmetry, and the two-parameter (g-deformcd) eight-state supersymmetric 
Q ' fermion model with quantum superalgebra Uq{gl{3\l)) symmetry, were introduced ||2l| , |2^ . 

O Dolcini and Montorsi [ p3[ introduced integrable extended Hubbard Hamiltonians from symmetric group solutions. 

^ One of the aims of this work is to show the solution of these models via the QISM. The most general form of the 
extended Hubbard model invariant under spin-flip and conserving the total number of electrons and magnetization, 
first considered in [e3, is described by the Hamiltonian 



ttST^ F v-^ t t 
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where fie is the chemical potential. Here, electrons on a lattice are described by canonical Ferimi operators Cj^a and 
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cj satisfying the anticommutation relations given by {c|^, Cj./j} = ^ijSap, where = 1, 2, and a, 13 ="1, ],. Also 

nj.a = Cj aCj^a, iT-j = JT-jj + (0) the term t represents the band energy of the electrons, while the subsequent 

terms describe their Coulomb interaction energy in a narrow band approximation: U parametrizes the on-site diagonal 
interaction, V the neighboring site charge interaction, X the bond-charge interaction, W the exchange term, and Y 
the pair-hopping term. Moreover, additional many-body coupling terms have been included in agreement with ]2^ : 
X correlates hopping with on-site occupation number, and P and Q describe three- and four-electron interactions. 

As was shown in |2^], the integrability of the model lies in the fact that there exists a solution of the Yang-Baxter 
equation which takes the form 



R{u) 



i + un 



where u is the spectral parameter and 
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Indeed, after making the assignment 
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|i)=|0), |2) =UT>,=44,tI0), |3) =IT),=4tI0)' |4> =U>,=4ilo> 

one can show that H can be expressed as the sum over a graded generalized permutator, 

L 



(2) 



H 



i=i 



(3) 



(4) 



where the operator flj^+i permutes the four possible configurations between the sites j and j + 1; viz. 

n = crfkieii Skk) + cTikieik <»s eki) 
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where crf/^ieu ekk) are diagonal terms and (7°k{eik <8) eki) are off-diagonal terms. At this 
mention that the construction here bears some resemblance to those by Alcaraz and Bariev | 
for solutions of the Yang-Baxter equation based on representations of the Hecke algebra. 

It is clear that this form of interaction conserves the individual numbers iVf and Ni of electrons with spin up and 
spin down respectively, and the numbers Ni and of doubly occupied (local electron pairs) and empty sites (holes) . 
We will choose the following conventions throughout this paper: 

A'^l=number of single electrons with spin up 

A'^l=number of single electrons with spin down 

Ne^Ni + A^i^number of single electrons 

A'i^number of local electron pairs 

A''?! = number of holes 

Nb=Nh + A^i=number of "bosons" . 

To be specific, we give the relations between tr^j,, cr°f. and parameters in the Hamiltonian (|l]): 
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af^=c, 42 = c+Q + U + 4P + 2fi,+W-2W, a^^ = af^ = c + V - W + 

d d , , U d d d d ,A^e 

cr^2 = (721 = C + + Cri3 = ^14 Cr31 = ^41 = C + y , 

4^ = 4a = <^32 = <^i2 = c + P +^ + ^ + 2V ~W, ai, = ai,^c + V + He 

It turns out that actually there are 96 different possible choices of values of the physical parameters in (1). They 
can be cast into six groups as (^) shows. 
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Above, there is the restriction si = ±1, i = 1, 5. 

We now construct the eigenstates of the Hamiltonians of the one-dimensional model in the above six groups, using 
the QISM. We use the R matrix satisfying the Yang-Baxter equation and introduce an L operator constructed directly 
from the R matrix of the twisted representation. The quantum Yang-Baxter equation can be written as the operator 
equation: 

R{X - /i)Lj(A) ® Lj{^) ^ Ljifi) ® Lj{X)R{X - fi). (5) 
Here (g) denotes the graded tensor product defined by 

where S Z2 denotes the grading of the index i. We chose to adopt the bosonic, bosonic, fermionic and fermionic 
(BBFF) grading ei = £2 = 0, £3 = 64 = 1 on the indices labelling the basis vectors. 

We now proceed to establish the relation between the Hamiltonian and the transfer matrix r(A), which is the 
supertrace of the monodromy matrix T(A) defined by 

T{X)^Ll{X)Ll-i---Li{X). 

From (|) it follows that 

R{X - n)T{X) (g) T{n) = T{p) (g) T{X)R{X - (6) 

Thus we have 

[T{X),T{fi)]=0 

and so the t(A) form a one-parameter family of commuting operators. The transfer matrices may be taken as integrals 
of the motion, and so we obtain an infinite number of higher conservation laws of the model. 

II. ALGEBRAIC BETHE ANSATZ FOR GROUP 1 

Having recalled the quantum integrability of the models, let us use the nested algebraic Bethe ansatz method to 
find the eigenvalues of the transfer matrices. As will be shown, the Bethe ansatz solutions for the six different groups 
take on differing forms. In particular, the number of levels of Bethe ansatz nestings ranges from to 2. Moreover, 
some cases do not admit a unique reference state. In such an instance, we are forced to use a subspace of reference 
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states to perform the calculations. This type of procedure was first investigated by Abad and Rfos 
of alternating sm(3) representations and we will adopt this formalism where necessary. 
We start from the first group with BBFF grading. The explicit form of the i?-matrix is 

R{u) = l + uU = 



for the case 
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To utilize the framework of the QISM, we write down the L operator 

1 



Lj{u) 



1 + SiU 



PR{u), 



where P is the graded permutation operator, 



(7) 



(8) 



We choose the local vacuum state as |0) j — (0, 0, 0, 1)*. Acting the L-operator on this local vacuum state, we have 





with a{u) = u/{l + siu). Define the vacuum state as |0) 
monodromy matrix in the following way: 
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i|0),. 



|0) 



(9) 



Using the standard QISM, we represent the 



T{u) = Ll{u)Ll-i{u) ■ ■ ■ Li{u) 



( A(u) Bi{u) B2{u) Bsiu) 

Ci{u) Dniu) Duiu) D^siu) 

C2{U) D2l{u) D22{U) D23{U) 

\C3{u) Dsiiu) D32{u) D33{u) 



The transfer matrix with periodic boundary conditions is thus written explicitly as 

t{u) = A{u) + Dn{u) - D22{u) - D^. 

The action of the monodromy matrix on the pseudo-vacuum state is 

Aiu)\0) = |0> Dniu)\0) = [53a(")]^|0) , D22iu)\0) = i?33(")|0) = [~a{u)]'^\0) , 
B,iu)\0)^0, Ck{u)\0) =0, Dkiium =0 (k^l, k, I ^1,2,3). 

Substituting ( ]lO| ) into the Yang-Baxter algebra (||) we may deduce the following commutation relations: 
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B,i\)Ai^i) - 4s.(A*)A(A)], 



a(A — n) 



Ba, {\)Ba, (/i) = r{\ - ^i)llZB,, ifi)B,, (A), 
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where 



r{u)^biu)lK: + aiu)uK: 



Here H*^^)"^ is the 9x9 submatrix of H, 
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where f i ' is a function of the spectral parameters Xj. Acting the transfer matrix (|11|) on the above vector, we 
have 

T{u)Ba, {\l)Ba, (A2) • • • Ba„ (A„)|0) 

= [A{u) + Dii{u) - D22{u) - D33iu)]BaA^i)Ba,{X2) ' ' ' Sa„(A„)|0) F'^"-'^^ 

n 

= Ss"' (-1)^= n (X _ ^ Ba. (A2) • • • Sa„ (A„)|0) F^^-'^^ 

j=i ^ i ' 

n ^ 

+ Ku)]^4-"(-l)^-+^^ n ^^ _ , ^'^n^,{Afc})g:;Xi3.,(Ai)i3,,(A2)---i3,JA„)|0)F°-----°^ 

j=i i 

+u.t. 

where u.t. denotes unwanted terms, and t'-^^(u) is the nested transfer martix. Denote the eigenvalues of t(u) and 
t(^)(u) by A(m) and A'^'^\u). We have 

In order to cancel the unwanted terms, we need the following Bethe ansatz equations 



The nested transfer matrix is written as the supertrace on the auxiliary space for the reduced monodromy matrix 
which satisfies the Yang-Baxter relation. 



t(i)(u, {Afc}) = str[diag(-, -1, -1)4'^ - A„)4'-i(^^ - A„_i) • • • l[^\u - Ai)] 

S3 

r(A - A^)r(i)(A) r(i)(/i) = r(i)(M) ® T(i)(A)r(A - f,). 



(13) 



If we write 
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TW(u)=iW(«)L«i(u)...L«(«) 

the L^^-' operator is 



/ b[^\u) B''i\u) 

I , 



where P'^^-' is the 9x9 permutation operator 

corresponding to the grading ei = 0, 63 = £3 = 1. Now ( ]l3| ) and r(M) imply that 

Bi\'> (A)P« (m) - r(i) (A - ^^t2<'' W- 

Here the values 1, 2 are both fcrmionic (ei = 1 = £2). The P-matrix r'^^^(/i) is 

r«(.):;=6(^)(.)/(^)::+a«(.)n(^):: 



1 + S5U 1 + S5U 



Above H^^^cd is a 4 X 4 submatrix of E^^^ 



n(2) = 



-S4 

-S4 

-S4 

-S4 



corresponding to the grading ei — €2 — I- As the reference state for the first nesting we choose |0) 
(1,0, 0)S |0) = (8)^^i|0) [^^ as the pseudo- vacuum. We find 

A^HIO) « = |0)«, 

n 

and due to t'^-'^'(m) = ^^'-'^^(m) + d[^^ (u) + g^t the eigenvalue 

jii 



A(i) («, { A,}) = - n — TT-Tn : + n — TT7 mT 11 ^^^^'^ - A^) A^^) (u, { A« }), 

provided the parameters {Am"*} satisfy 



^3 fJi S2a(i) (A W - u) fJi S2a(i) ~ A W ) ' 



Z = l 

The transfer matrix of the second nesting is written as 



{AW}) = .tr[L(2)(^ - AW)Piti(- - Ai?-i) • • • if - A^)], 



where 

~h'^^\u)ef a(2)-&(2)(u)ef j • ^^^^ 
From the Yang-Baxter relation for t^^^(m) one can derive the foUowing intertwining relation 

r(i)(A - A*)r(?)(A) ® Tif{^) = r(2)(^) ® T(f (A)r(i)(A - /i). (18) 



The components of (18) needed for the construction of an algebraic Bethe ansatz are 



B^^\\)B^^\^i) = i?(2)(^)ij(2)(A), (19) 



where 



1 



1 + SiU 1 + SiU 

For the reference state for the second nesting we pick |0) = (1,0)*, |0) — ''^^ ■ From the action of the 

nested monodromy matrix 

Tif(u)=L(^)(.)<u(.)...Lf)(.). 



we find 



j=i -^(Aj - u) 



due to r(2)(u) = -^(2)(u) -D(2)(u). Thus 



A^^'(.,{A«}) = -in ^iyT^ ft + n , n w^^.- a^)) 

,1 540^-^^ (w -K )m=i o-^ ' (Am - m) ,=i S4a^2) (A ■ - it) 



under the condition that the spectral parameters {Ap^''} are solutions to the Bethe ansatz equation 

n \ u \ = n ^4«^''(A« - Af ). 

We have now obtained the complete set of nested Bethe ansatz equations, which read 
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11 ^ ,>(2) ,(2)^ , -, ^ 11 
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Here we have used n = + Ni = + Ni + Ni,ni = Ne,n2 = Ni. The corresponding energy eigenvalue E of the 
model is given by 



1 



L. 



(20) 



The BBFF grading solution of above models closely follows the solution of the EKS model in |p|-|ri|]. 



III. ALGEBRAIC BETHE ANSATZ FOR GROUP 2 



The algebraic Bethe ansatz calculations for this group proceed in exactly the same manner as group 1 up to the 
introduction of the matrix II'^'. which for this case reads 
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Following the calculation along the same lines as the previous section, we find that the matrix (u) appearing in 
the equations is of the form 

Here /^^■'aa is the 4x4 identity matrix. For the reference state of the first nesting we choose the state |0) ^f^'' — 
(1, 0, 0)*, |0) = ®fe=i|0) as the pseudo-vacuum and then find 

n 



Due to t(i)(u) = + d[\\u) + D^^^{u), we get the eigenvalue 
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n 



under the condition that the spectral parameters {Am''} are solutions of the Bethe ansatz equation 

fc=l 
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For the full solution, we have the nested Bethe ansatz equations 
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where n, ni, 77.2 have the same meaning as previously. Also, the energy expression (]20|) applies here 



IV. ALGEBRAIC BETHE ANSATZ FOR GROUP 3 



We now consider the case of the algebraic Bethe ansatz for group 3. As we will see here, the procedure is funda- 
mentally different from the preceding cases in that we are required to work with a subspace of reference states for the 
first level of the algebraic Bethe ansatz. The methodology we employ follows that proposed by Abad and Rios pTj . 

In the case of group 3, the i?-matrix reads 
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and we again express the L operator as 



If we choose the local vacuum state as |0) j 
we have 
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:(a, /3, 0, 0)*, and act the L-operator on this local vacuum state, 
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and represent the monodromy matrix as 



T{u) = Ll{u)Ll-i(u) ■ ■ ■ Li{u) 



/ Aii{u) Ai2{u) Bii{u) Bi2{u) 

A2l{u) A22(m) B2i{u) B22{u) 

Cii{u) Ci2(u) Dn(u) D^2{u) 

\C2l{u) C22N D2l{u) D22{u) 



(23) 



The transfer matrix is thus written explicitly as 

t{u) = Aii{u) + A22{u) - Dii{u) - i:'22(w)- 

The action of the monodromy matrix on the vacuum state is 

[Aii{u) + A22(m)]|0) = tro[PwPL-ifl ■ ■ ■ Pio]|0) 

i?ii(«)|0) ^D22{um ^[a{u)]\-ir^s^\ 
B,fc(w)|0)^0, Gfc(u)|0)=0 (^,A: = 1,2), 

where Pjq is the permutation operator for two-dimensional spaces (corresponding to the indices 1 and 2). Substituting 
(p3|) into the Yang-Baxter algebra (||) , we may deduce the following commutation relations: 
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a(A — /ij a(A — /ij 

5ac(A)BM(M) = '^(A - ll)i'^Bac'{f^)BM'W, (24) 

where 

Here H^^^"^ = ssF^^^ with permutation matrix P^^^ = — ^ Cy (g) Cji corresponding to the grading ci = £2 = 1 and 
Si = —1, 5*2 = j^- Denote the eigenvalues of t(u) and t'^^^ by A{u) and A(^)(m). We now have 

i r / N T r TT i 



Here G = tr[diag(— 1, j^)PloPl-i.q ■ ■ ■ Pio], and the parameters {Afe} are required to satisfy the Bethe ansatz equations 

ni 



A(i)(A,) = [a(Afe)](-^).G-[] 



a(Afc - Aj) 
_^ «(Aj - Afc) 

The nested transfer matrix is written as the supertrace on the auxiliary space for the reduced monodromy matrix 
which satisfies the Yang-Baxter relation, i.e. 

TW(u,{Afc}) = str[diag(-l, -)Li'^ {u - X„,)Li%{u ~ X^.^i) ■ ■ ■ l['\u - X^)], 

r(A - m)T(^)(A) ® T(i)(m) = r,W(M) ® T(i)(A)r(A - ^,). (25) 



The components of (25) needed for the construction of an algebraic Bethe ansatz are 



B(i)(A)B(i)(/i) =B(i)(^)b(i)(A), 



where 
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As the reference state for the second nesting we take |0) = (1,0)*, |0) = <8)feii|0) i^'. From the action of the 
nested monodromy matrix 



we find 
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and due to t^'^^u) = A'-'^^u) - —D'-^^u) we have 
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M^3a(^H^^-A«)H«^'HA.-H) . li^^^(3)(Ai)_,)li/ 
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under the condition that the spectral parameters {A„i } are solutions to the Bethe ansatz equation 

Ua(3)(A«-Af))- ^^M' 
Now we obtain the complete set of nested Bethe ansatz equations reading 

/ l + giAfc \'-_ 1 ^ ,si(Afc - A,) - 1 gi(Afc - A;) - 1 g3(Af ^ - \k) - 1 

i^ ,s3(Ay^-AW)-l _ if^ .3(A.-A<^)) 
The energy expression for this model reads the same as in the previous cases (|20|). 



V. ALGEBRAIC BETHE ANSATZ FOR GROUP 4 



Just as the calculations for the cases of group 1 and 2 follow along similar lines, we find an analogous situation 
occurring with groups 3 and 4. For group 4, we have the i?-matrix 



R{u) 



ull 



+ SlU 















































1 


+ Silt 


















































1 

















~u 
































1 


























—u 




















1 


+ Silt 


















































1 + Silt 


















































1 








S2U 









































1 

















S2U 































1 









































S21t 








1 


















































1 + Salt 


















































1 + S3?t 























~u 


























1 

















































1 















































1 


+ S3li 


















































1 + S3li 



(26) 



V 

The calculations of the algebraic Bethe ansatz proceed in exactly the same manner as the group 3 case except now 
we find that in (pih we have 



bb 



1 + S3U 

r 



/(2) 



bb 



SiU 



For the eigenvalues A(u) of r(u) we obtain the expression 

rii Ml 



a{u-\j)'' 



a(Xj — u) 

where G is defined as before and {A^} are subject to the Bethe ansatz equations 

'l + si\k\^ ^ -PT ^i(Afc - Aj) + 1 



Afc 



G n 



si(Afc - Xj) - 1' 



Again, the energies are given by (EG). 
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VI. ALGEBRAIC BETHE ANSATZ FOR GROUP 5 



For group 5, the i?-matrix is given by 
R(u) = 1 + tin = 



+ Sill 


















































1 








S2U 









































1 

















— tt 
































1 


























— tt 














S211 








1 


















































1 + S'iU 


















































1 + Salt 















































1 


+ S3tt 
































—u 

















1 















































1 


+ Sstt 


















































1 + S3tt 















































1 


+ S[iU 























—u 


























1 















































1 


+ SgU 


















































1 + S3tt 


















































1 + satt 



(27) 



v 

We choose the local vacuum state as |0) — (1, 0, 0, 0)*. Acting the L-operator on this local vacuum state, we have 



1 











S2a{u) 




























-a(ti 



|0) 



(28) 



Define the vacuum state as |0) — <8'5'=i|0) -. The monodromy matrix is represented as 



T{u) = Ll{u)Ll-i{u) ■ ■ ■ Li{u) 



A{u) 
Ci{u) 
C2{u) 



Bi{u) 
Dn{u) 
D2i{u) 



B2{u) 

Dl2(u) 



B2{u) 



£'22 (u) D23{u) 
Dz2{u) D33{u) 



(29) 



and so the transfer matrix is explicitly 

t{u) = Aiu) + Dn{u) - D22[u) - Dsaiu). 

The action of the monodromy matrix on the vacuum state is 

v4(tt)|0) - |0), Z?ii(ii) = [S2a(w)]^|0), D22iu) = D33iu) = [-aiu)]'^\0) 
Bkiu)\0)^0, Ckiu) = 0, Afe(u) = 0, {i^k, i,k=l,2,3). 

Substituting ( p9| ) into the Yang-Baxter algebra (||) , we may deduce the following commutation relations: 



(30) 



Dab{^i)B,{X) = Sa[{-1) 

^(m)Bc(A) 



Ba.WBaAt^) 



a{\ — /i) 
r(A-At)^?^?B, 



a{fi — A) 
B,(A)A(/i) 



''''' Bt{x)DM - {-ly 



a(/i — A) 



Bb{fl)Dac{X)], 



a{X — /i) 



Bci^^)AiX)], 



ai ai -^ai 



where 



bb 



1 + S3U 
1 + SiU 



/(I) 



bb 



Here Z^^-'at the 9 x 9 id matrix and 5*1 = -^,82 = —1, S3 = —1. The eigenvalues of t(u) read 



A{u) 



-N, 



i-ir- n 



L a{Xj - u) 



Hu)] 



n 



\ a{u- X^) 
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with the following Bethe ansatz equations 

L 



1 + siA, 



-We „L-n 



(-1) 



n 

fc=i 



si{Xj - Afe) ~ 1 



and the energy is given by (E 



VII. ALGEBRAIC BETHE ANSATZ FOR GROUP 6 



The final case to consider corresponds to the i?-matrix 
R{u) = 1 + un = 



1 + SlU 















































\ 





1 








S2U 









































1 + Silt 

















































1 


-f Silt 











































S2U 








1 




















































1 + Slit 




















































1 








u 











































1 

















u 































1 


-f Silt 











































u 








1 




















































1 + Slit 




















































1 + SlM 

























—u 























1 




































u 

















1 

















































1 


-f Slit 




















































1 + SlM 


/ 



V 

In contrast to the other cases considered, we choose the local vacuum state as |0) j 
on this local vacuum state, we have 



(31) 



(0,1,0,0)*. Acting the L-operator 



S2a{u) 

* 1 * 

a{u) 






(32) 



Defining the vacuum state as |0) 



®f=i|0) express the monodromy matrix as 



T{u) = Ll{u)Ll~i{u) ■ ■ ■ Li(u) 



/ Dn{u) Ci{u) Di2{u) Duiu) 

Bi{u) A{u) B2{u) B2{u) 

D2l{u) C2{u) D22{U) D23iu) 

\D3i{u) Csiu) D32{u) Dssiu) 



(33) 



and so the transfer matrix is 



t{u) = Dn{u) + A{u) - D22{u) - D^^iu). 

The action of the monodromy matrix on the vacuum state is given by 

Diiiu) = [S2a(u)]^|0) , A{u)\0) = |0) , D22{u) = D^siu) - [a(w)]^|0) 
Bk{u)\0) ^0, Cfc(u) = 0, Afe(w) = 0, (^^fc, i,A:=l,2,3). 

Substituting (|33| ) into the Yang-Baxter algebra (||) we find 



(34) 



a(3)(/i- A)' 



a(3)(^- A)" 



Ai^l)BciX) - Se 



1 



V(3)(A-m) a(3)(A-/i) 
BMBM = r{X - y)lTalBaMBaAX), 
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where 

i + S3U 

and now S'l = 52 = 1, = 1. The eigenvalues for the transfer matrix read 

^(") - -2"''' n „(3)(A^.„„) + [«(")]^-2 n ai^Ku-X,) 
so that the Bethe ansatz equations 

^2 «2 11 



fc=i 



are satisfied. Here the energy eigenvalue differs somewhat from the previous cases and has the form 

L-Ni 



^ A,(l + siA,) 



VIII. SUMMARY AND DISCUSSION 



In this paper, integrable extensions of the Hubbard model arising from supersymmetric group solutions, by means 
of the algebraic Bethe ansatz method, have been investigated. In particular, we have calculated explicitly the Bethe 
ansatz equations as well as the energy eigenvalues for 6 different classes of underlying i?-matrices, which in fact 
correspond to 96 different possible physical Hamiltonians. 

A natural direction for possible further research is to deal with physical applications of the above models. More 
specific future works will be: (i) studying low energy behaviour and physical properties of the corresponding systems 
based on an analysis of the Bethe ansatz equations from these results, including: investigating the ground state 
structure, computing the finite size corrections to the low-lying energies, and calculating thermodynamic equilibrium 
properties, using the methods of Woynarovich (ii) employing some traditional mathematical methods such as the 
Wiener-Hopf technique to solve the special kind of integral equations arising from the thermodynamic Bethe ansatz 
equations, using the methods of Yang & Yang and Babujian p9| . 



This work is supported by the Australian Research Council. 
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